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Using molecular dynamics simulations, we study dynamics of a model polymer melt composed 
of short chains with bead number N = 10 in supercooled states. In quiescent conditions, the stress 
relaxation function G(t) is calculated, which exhibits a stretched exponential relaxation on the time 
scale of the a relaxation time r a and ultimately follows the Rouse dynamics characterized by the 
time tr ~ N 2 r a . After application of shear 7, transient stress growth o- xy (t)/j first obeys the linear 
growth j*dt'G(t') for strain less than 0.1 but saturates into a non-Newtonian viscosity for larger 
strain. In steady states, shear-thinning and elongation of chains into ellipsoidal shapes take place 
for shear 7 larger than r^ 1 . In such strong shear, we find that the chains undergo random tumbling 
motion taking stretched and compact shapes alternatively. We examine the validity of the stress- 
optical relation between the anisotropic parts of the stress tensor and the dielectric tensor, which 
are violated in transient states due to the presence of a large glassy component of the stress. We 
furthermore introduce time-correlation functions in shear to calculate the shear-dependent relaxation 
times, T a (T, 7) and Tr(T, 7), which decrease nonlinearly as functions of 7 in the shear-thinning 
regime. 

PACS numbers: 83.10.Nn, 83.20.Jp, 83. 50. By, 64.70.Pf 



I. INTRODUCTION 
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The dynamics and rheology of glassy polymers are 
known to be very complicated and are still not well un- 
derstood. We summarize salient features of such systems 
below. 

First, in the linear response regime, thermal relax- 
ations of the chain conformations occur from microscopic 
to macroscopic time scales, as repealed in measurements 
of stress and dielectric responses.Ero In a relatively early 
stage, the stress relaxation function G(t), which describes 
linear response to small shear deformations, can be fitted 
to the Kohlrausch- Williams- Watts (KWW) form, 



G G (i)=G exp[-(i/T s ) c ], 



(1.1) 



after a microscopic transient time it ra - The time t s (3> 
T tra ) is of the order of the structural a relaxation time 
r a (to be defined in (3.8) below), which grows dramati- 
cally as the temperature T is lowered towards the glass 
transition temperature T s . When the time t consider- 
ably exceeds r s , the relaxation of the chain conforma- 
tions is relevant and is well described by the Rouse or 
reptation dvjpamics, depending on N < N e or N > N e , 
respectively.u Here N is the polymerization index and 
N c is that between entanglements on a chain. For short 
chain systems with N < N c , the overall behavior in the 
time region t 3> ftra may be expressed as 



G(t) = G G (t) + G R (t). 



(1.2) 



The C?R,(t) is the stress time-correlation function in the 
Rouse model whose terminal relaxation time tr is of or- 
der N 2 r a . For entangled chain systems with N > N e , 
the KWW function in (1.1) is followed by the power-law 
decay, 



with v rsj 0.5 until the rubbery plateau G(t) 



reached,Bi where G^ 1 assumes the modulus nk-QT/N c 
of entangled polymers with n being the bead number 
density. Ultimately, G(t) follows the reptation relax- 
ation G lep (t) on the time scale of a very long rep- 
tation time T rcp . These hierarchical relaxations arise 
from rearrangements of jammed atomic configurations 
and subsequent evolution of chain conformations. They 
also give rise to the corresponding characteristic behav- 
iors in the frequency-dependent shear modulus G*(u>) = 
ioj / °° dte~ lu)t G{t), depending on the frequency uo relative 
to the inverse characteristic times introduced so far.B~H 

Second, in the nonlinear response regime, glassy flu- 
ids generally exhibit highly viscous non-Newtonian flow 
close to (but .above) T g even if they are low-molecular- 
weight fluids. □ In such fluids, if 7 > t^ 1 , atomic rear- 
rangements are induced not-.tw thermal agitations but 
by externally applied shear.ETLl In chain systems with- 
out entanglements, on one hand, shear thinning occurs 
at sufficiently high (but sometimes uraeaiistically large) 
shear rates due to chain elongation.BlI3 In entangled 
polymers, on the other hand, shear thinning occurs at 
very small shear largerthan r~i , where disentanglements 
are induced by shear.H Thus supercooled chain systems 
are most easily driven into a nonlinear response regime 
even by extremely small shear, though the crossover shear 
stress from linear to nonlinear regimes may not be very 
small. Furthermore, in glassy fluids below T g , plastic de- 
formations are often induced in the form of large-scale 
shear bands, .above a yield stress (corresponding to a few 
% strain) Jl3'EJ It is of great importance to understand 
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how these nonlinear effects occur depending on 7, T, and 
N. 

Third, in rheological experiments on polymers, use has 
been made of the stress-optical relation between the devi- 
atoric (anisotropic) parts of the dielectric tensor e Q/ 3_(ai 
optical frequencies) and the average stress tensor er a/ 3£rt£l 
In shear flow with mean velocity jy in the x direction, it 
is expressed as 



Cq(<t 



crpp) 



(1.4) 



where Cq is called the stress-optical coefficient. For melts 
this relation excellently holds at relatively high T (> T g ) 
for general time-dependent nonlinear shear deformations. 
If measurements are made in steady states, it holds even 
close to Tg.EEl For its validity, we need to require that 
the form contribution to is negligible as compared to 
the intrinsic contribution!!.! and that the glassy part of 
the stress is negligible as compared to the usual entropic 
part. Thus it is violated when the form part is relevant 
such as in polymer solutions close to the demixing critical 
point or when measurements are made in transient states 
close to T g . In the latter case, as is evident from the en- 
hancement of t s in (1.1), the glassy part of-ihe. stress is 
dominant for relatively rapid deformations. t3~ll3 

While the predictive power of analytic theories_in_poly- 
mer science is still poor, computer simulationscj~E3 can 
provide us a useful tool to investigate the microscopic 
origins of experimentally observed macroscopic phenom- 
ena. In quiescent states, diffusive motions in supercooled 
melts have beeruextensively studied-using molecular dy- 
namics (MD)ErE3 and Monte CarloE3~El simulations. In 
another application, nonequilibrium molecular dynam- 
ics (NEMD) simulations have been usefjal[4flj»Mestigate 
chain deformations and rheolog y in flow.^BBEl In par- 
ticular, Kroger et alEU studied the molecular mechanisms 
of the violations of stress-optical behavior for a melt con- 
sisting M — 260 chains with bead number N — 30 after 
application of clongational flow. 

In this paper, we will present results of very long MD 
simulations to study linear and nonlinear dynamics of a 
supercooled polymer melt in the absence and presence of 
shear flow. Long simulation times are needed to calculate 
the terminal relaxation of G(i), which has not yet been 
undertaken in the literature. As a new finding, we will 
show that each chain in our melt system is changing its 
orientation (tumbling) randomly in shear flow. Use will 
be made of techniques and concepts introduced in our 
previous papers on supercooled binary mixtures under 
shear flowxJD Some of our results were already published 
elsewhere. Qc3 



II. MODEL AND SIMULATION METHOD 

Our system is composed of M — 100 chains with 
N = 10 beads confined in a cubic box with length 
L = 10cr and volume V = L 3 = 10 3 er 3 . The number 



density is fixed at n = NM/V = 1/er 3 , which results in 
severely jammed configurations at low T. All the bead 
particles interact with a truncated Lennard- Jones poten- 
tial defined byE3 



U hJ (r) =4e 



+ e {r<2 1 ' & a). (2.1) 



The right hand side is minimum at r = 2 1 / 6 er and 
the potential is truncated for larger r (Uz,j(r) = for 
r > 2 1/6 er). By using the repulsive part of the Lennard- 
Jones potential only, we may prevent spatial overlap of 
the particlesHl Consecutive beads on each chain are con- 
nected by an anharmonic spring of the form, 



Mr) = ~\k c Rl\n[l - (r/R ) 2 } 



(2.2) 



with k c = 30e/er 2 and Rq = 1.5a. In our simulation the 
bond lengths b) = \R) - R% +1 \ (1 < j < N - 1) be- 
tween consecutive beads on the same chain k were very 
close to the minimum distance 6 m i n = 0.96ct of the sum 
Uh.i(r) + Uf{t). The deviations 6^ — b m i n were only on 
the order of a few % of 6 m i n for any T and 7 realized in 
our study. 

Microscopic expressions for physical quantities can be 
expressed in terms of the momentum and position vec- 
tors of the j-th bead on the k-th chain, and p 1 *, where 
j = 1, • • • , N and k = 1, • • • , M. For example, the space 
integral of the microscopic stress tensor reads 

all pairs 



k=lj=l 
M N-l 



k=i j=i 



(2.3) 



where m is the mass of a bead, U^AO = dUi,j(£)/cl£, 
and £7^(0 = dU F ($/dt. Here £ = &) in the right 



hand side represents the relative vector _R 7 —R4, between 



the two beads, Rj and Rj, , in the second term and the 
relative vector R^ — Rj + i between the two consecutive 
beads, Rj and Rj+i, of the same chain in the third term. 
To avoid cumbersome notation, we will write the bead 
positions simply as Rj (j — 1 , ■ ■ • , N) suppressing the in- 
dex k. When they will appear in the statistical averages 
(•••), the average over all the chains Ylh=i(' ' )/^ wu ^ 
be implied even if not written explicitly. Furthermore, 
it is convenient here to introduce the usual notation o a p 
for the stress tensor by 



1 T 



J af3, 



(2.4) 



where p is the pressure and the second term is deviatoric.i 
The a a {} has already appeared in the stress-optical rela- 
tion (1.4). 
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We will hereafter measure space and time in units of 
a and r = (mcr 2 /e) 1 / 2 . The temperature T will be mea- 
sured in units of e/fce ■ The original units will also be used 
when confusion may occur. Our simulations cover nor- 
mal (T — 1.0) and supercooled (T = 0.2) states with and 
without shear flow (7 = 0, 10" 4 , 10~ 3 , 10~ 2 , and lO" 1 ). 
Simulation data were taken after very long equilibration 
periods (~ 10 2 t r ~ 5 x 10 6 at T = 0.2) so that no ap- 
preciable aging (slow equilibration) effect was detected in 
the course of taking data in various quantities such as the 
pressure or the density time-correlation functions, (i) In 
quiescent cases, we impose the micro-canonical condition 
and integrate the Newton's equations of motion, 



dt 1 



1 J 



dt P > 



fi 



(2.5) 



where fj is the force acting on the particle j due to the 
potentials. Integration was performed with time incre- 
ment At = 0.005 under the periodic boundary condition. 
Long time simulations of order 10 2 tr, which corresponds 
to 10 9 MD steps fot-X-j=r Q.2. pere performed. In the pre- 
vious simulationsJ^Sl^EaLy however, the integrated 
times did not much exceed tr in supercooled states, (ii) 
In the presence of shear, rewriting the momentum devi- 
ations Pj — m A jYje. x from the mean flow as Pj, we, .inte- 
grated the so-called SLLOD equations of motion,! 



d - 



(2.6) 



where e x is the unit vectors in the x (flow) direction, 



(Xj,Yj,Zj), and p = (p x j,Pyj,Pzj)- 



coefficient ( was set equal to 



The friction 



(2.7) 



The temperature T (= ^ ■ p^/mNM) could then be kept 
at a desired value. The time increment was At = 0.0025. 
After an equilibration run in a quiescent state for t < 0, 
we gave all the particles the average flow velocity jYje x 
at t = rajld, then imposed the Lee- Edwards boundary 
conditionEj'E 4 ] to maintain the shear flow. Steady sheared 
states were realized after transient relaxations. 



T p = (b 2 /[l2k B Tsm 2 {ir P /2N)}, 



(3.1) 



where 1 < p < N — 1. The Rouse relaxation time tr is 
the slowest relaxation time, 



tr = n N 2 (b 2 /(3ir 2 k B T). 



(3.2) 



The segment length b in the corresponding Rouse model 
may be related to the variance of the end-to-end vector 
of a chain P = Rn — Ri in our microscopic model by 



\P\ 2 ) = b 2 {N-l). 



(3.3) 



As a result, b is dependent on T but its dependence turns 
out to be weak as b = 1.17, 1.18, 1.19 for T = 1.0, 0.4, 
0.2, respectively. Note that b is larger than the mini- 
mum distance 6 m i n = 0.96 of the bond-potential. Let us 
consider the time-correlation function of P(t), 



C(t) = (P(t + t )-P(t ))/(\P\ 2 ), 



(3.4) 



which is normalized such that C(0) = 1. Here C(t) 
should be independent of the initial time to in steady 
states in the limit of large system size. However, our 
system is not very large, so we took the average over 
the initial time to. This statistical averaging will not be 
mentioned hereafter in showing our MD results of time- 
correlation functions. In the Rouse dynamics C{t) is cal- 
culated as 



C R (t) 



(N - l)N ^ 

v ' odd p 



cot 2 



irp 
2N 



-t/Tp 



(3.5) 



where the summation is over odd p but the first term 
(p = 1) is dominant in the whole time region (so we may 
determine tr by C(tr) = e" 1 ). As shown in Fig.l, our 
MD data of C(t) can be excellently fitted to C R (t). The 
tr thus determined increases drastically with lowering T 
as tr = 250, 1800, and 6 x 10 4 for T = 1.0, 0.4, and 0.2, 
respectively. In.the.preinons simulations on nonentangled 
polymer melts, E 4 Ji 2 J : E2l~OE3 numerical results were con- 
sistent with the Rouse dynamics for small p (large-scale 
motions), but deviations are enhanced for large p (small- 
scale motions) in supercooled states. Furthermore, we 
here give the expression for the stress relaxation function 
in the Rouse model, 



III. DYNAMICS IN QUIESCENT STATES 



s been confirmed 
that the single- 



Although it is highly 
by computer simulationa 2 . 1 
chain near-equilibrium dynamics in unentangled melts 
can be reasonably well described by (or mapped onto) 
the simple Rouse model. In the Rouse dynamics, the re- 
laxation time of the p-th mode of a chain is expressed in 
temis of a friction coefficient £ and a segment length b 

aska 



G R (t) 



nk B T 
N 



N-l 

E 

P =i 



exp(-2t/T p ), 



(3.6) 



which is equal to nk B T(N — 1)/N at t = and de- 
cays as nksTN -1 exp(— 2t/r^) for t > tr. Since G(t) 
is much larger than Gr(£) in the relatively short time 
region t < t Si they can coincide only in the late stage. 

In Fig||, we show the van Hove self-correlation func- 
tion, 
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^(^^Hig-Afl^f, (3.7) 

3=1 

where q = 2ir, ARj(t) = Rj(t + to) — Rj(to) is the dis- 
placement vector of the j-th bead in the time interval 
[to, to + 1] . The peak wave number of the static structure 
factor is given by q = 2ir. We define the a relaxation 
time r a from the condition, 

FM^e" 1 (<7 = 27r). (3.8) 

As has been reported in thejiterature, r Q increases dras- 
tically with lowering T.ExBHEl In our case, we ob- 
tained r a = 0.91, 5.8, and 310 for T = 1.0, 0.4, and 
0.2, respectively. At T = 0.2, where the particle motions 
are considerably jammed, F q (t) exhibits a two-step re- 
laxation and may be excellently fitted to the stretched 
exponential decay (oc exp[— (t/-r Q ) - 64 ]) for t > 10. Thus 
our system at T = 0.2 has characteristic features of a 
supercooled state, although its melting temperature is 
unknown. We find r Q ~ 10~ 2 (b 2 /k R T and tr ~ N 2 r a at 
any T. Particularly, for T = 0.2, we obtain 

T a =i 0. 017 C,b 2 /k B T, tr =i l.9N 2 r a . (3.9) 

The friction coefficient ( in the mapped Rouse model 
grows strongly as T is lowered in supercooled states. 

Now we discuss the linear viscoelastic behavior in su- 
percooled states. In terms of U xy (t)Jm(2.3), the stress 
relaxation function G(t) is written ascjfiJ 

G(t) = (Uj y (t + t )I$ y (t ))/k B TV. (3.10) 

In Fig.|[ we show numerical data of G(t), where the av- 
erage over the initial time to was taken but the data 
become noisy at very large t ~ tr. In the very early 
stage t < 1, G{t) oscillates rapidly due to the vibra- 
tions of the bond vectors bj = Rj — Rj+i. The initial 
value G(0) takes a large value (~ 100 in units of e/a 3 ) 
nearly independent of T. For T = 0.2, G(t) can be nicely 
fitted to the the stretched exponential form (1.1) with 
Go = 5, t s = 90 = 0.33r Q , and c = 0.5 in the time region 
1 < t < 10r s . For t > 50r s it approaches the Rouse stress 
relaxation function Gr(£) in (3.5). The zero-frequency 
Newtonian viscosity is given by 77(0) = / °° dtG(t), so it 
consists of the glassy (monomeric) part, 

poo 

m = / dtG Q (t) ~ 10r s , (3.11) 
Jo 

and the Rouse (polymeric) part, 

/>oo 

r /R = / dtG R (t) S O.SOSTAT-Vr. (3.12) 
Jo 

The ratio jyc/^R is of order 1/TN. They are of the same 
order in the present case of N = 10 and T = 0.2. How- 
ever, we should have rye <C 7/r for much larger N. 



To examine the orientation of the bonds, we consider 
the orientational tensor, 

chain j — 1 

where b^bj are the normalized bond vectors since |6j| = 
bo as stated below (2.2). Notice that in the Rouse model 
the space integral of the entropic stress tensor is given 
by the expression o\p = {ik^,TbQ/b 2 )Q a p, where b de- 
termined by (3.3) appears instead of bo- To compare our 
microscopic system and the simplified Rouse model, we 
calculated the time-correlation function, 

G h (t) = (a h xy (t + to)a h xy (t ))/k B TV, (3.14) 

by integrating (2.5). As shown in Fig.^J, Gh{t) is fairly 
close to Gr(£) in (3.5) from the Rouse model. In partic- 
ular, for t > O.Itr, we find G(t) =i G b (t) =i G R (t). 



IV. STEADY STATE BEHAVIOR IN SHEAR 
FLOW 

In Fig.[|, we display the steady-state viscosity 77(7) = 
&xy/i obtained at T — 0.2, 0.4, and 1, where the time 
average of the stress was taken. The crossover shear rate 
from Newtonian to shear-thinning behavior is given by 
T R 1 ~ N^ 2 ^ 1 . We may introduce the Weisenberg num- 
ber Wi by 

Wi = 7tr - 77V 2 r Q . (4.1) 

In the non-Newtonian regime, we have Wi > 1. The 
shear stress at the crossover is of order nk^TN^ 1 , which 
is the elastic modulus of the Rouse model. The hori- 
zontal arrows indicate the linear Rouse viscosity t?r in 
(3.12), while the vertical arrows indicate the points at 
which 7 = Tj^ 1 . In particular, the curve of T = 0.2 may 
be fitted to 

77(7) oc r v (4-2) 

with v = 0.7 for 7TR > 1. The 77(7) becomes insensitive 
to T for very high shear (Wi 3> 1). In MD simulations 
of short chain systems in normal liquid states,cTtL3 simi- 
lar shear-thinning has been reported, where the crossover 
shear is much higher, however. InJVLD simulations of su- 
percooled simple binary mixtures J3~l3 shear-thinning be- 
comes apparent for 7 > t^ 1 , where r Q grows at low T 
and which is consistent with (4.1) if we set N = 1. 

To demonstrate the stress-optical law in steady states, 
we show steady-state data of a xy /T vs Q xy for T = 
0.2,0.4, and 1 in Fig. 6. If the electric polarization ten- 
sor of a bead is uniaxial along the bond direction, the 
deviatoric part of the dielectric tensor is proportional to 
that of ther-tensor Q a p in (3.13). In accord with the 
experimentJiS our data collapse onto a universal curve 
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independent of T both in the linear (Q xy < 0.05) and 
nonlinear (Q xy > 0.05) regimes. 

We next consider anisotropy of chain conformations 
in shear flow. In Fig.||(a), we plot the x-y cross-section 
(z = 0) of the steady state bead distribution function, 

i N 

gs(r) = jjY.( S ( R i- R v-r))> ( 4 -3) 

3=1 

where 7 = 10~ 4 , T = 0.2, and R G = N' 1 ~£f=i Rj is 
the center of mass of a chain. In Fig.|?](b), we also plot 
the structure factor in the q x -q y plane (q z = 0), 

1 N 

which is proportional ±0 the scattering intensity from la- 
beled chains in shear. l3 In these figures 9 is the relative 
angle of the ellipses with respect to the y (shear gradient) 
direction. These figures demonstrate high elongation of 
the chains for 7 > t^ 1 . As will be shown in Fig.|^ below, 
they almost saturate into the shapes shown in Fig.|?] once 
7 exceeds r^ 1 . 

Let us define the tensor, 

j JV N 

^ = Jp - Rjc)(Ri0 - Rip)) 

»=i 3=1 

2 N 

= n - R G a )(RjP - Rgp)). (4.5) 

3=1 

For small q = (q x ,q y ,0), S(q) is expanded as 
S(q) = l-~ E / «/39a?/3H 

= l-i a ?(q- ei ) 2 -ia2( 9 .e 2 ) 2 + ..., (4.6) 

where {ei, and {a 2 , a 2 .} are the unit eigenvectors and 
eigenvalues of the tensor I a p (a, (3 € £,y). The two 
lengths ai and et 2 correspond to the shorter and longer 
radii in the principal axes of the ellipses. In terms of 
6, we have e\ = (— sin 6*, cos 9) and e 2 = (cos 9, sin 9) in 
the x-y plane. In Fig|], we display tan 9 = —e.\ y je\ x , 
1 — 01/02, and the xy component of the alignment tensor 
Q xy in (3.13). All these quantities represent the degree of 
deformations of chain conformations in shear flow. They 
are insensitive to T if plotted vs 7TR,. For 7TR, < 1, tan# 
is close to 1 (9 = 45°) and both 1 — ai/a 2 and Q xy lin- 
early increase with increasing 7TR. For 7TR, ^ 1, these 
quantities saturate into limiting values. At T = 0.2, they 
are 

(9^80°, oi/a2 = 0.3, Q xy = 0.1. (4.7) 
These are consistent with Q xy ~ sin cos 0. 



V. TRANSIENT VISCOELASTIC BEHAVIOR 

In Fig.^J, we plot the viscosity growth function a xy (t) /*/ 
after application of shear 7 at t — for T = 0.2. The 
system was at rest for t < 0. In the initial stage jt < 0.1, 
it evolves following the linear viscoelastic growth, 

1 r* 

-a xy (t) = / dt'G(t'). (5.1) 
7 Jo 

In the nonlinear regime, <J xy {t)/^ tends to the non- 
Newtonian viscosity 77(7). As a guide, we also display the 
linear growth function J Q dt'Gn(t') in the Rouse model. 
In the very early time region 1 « ( < r a , the growth 
(= Got) is much larger than the Rouse initial growth 
(= kftTt). We can see a rounded peak for 7 = 0.1 be- 
fore approach to the steady state in Fig.|[ More pro- 
nounced overshoot was already reported at high shear in 
MD simulation of much longer unentangled alkane chains 

(Cl{X)H202).li3 |— - |— - 

In experimentsJiTliia the stress-optical relation is tran- 
siently violated at low T after application of elongational 
flow due to the enhancement of the glassy component 
of the stress. For shear flow, Fig.|o| displays our MD 
results at T = 0.2 after application of shear at t = 
in a stress-optical diagram, where the solid lines are the 
averages over ten independent runs. As time goes on, 
the system traces the curve of a given 7, passes across 
the dashed curve representing the steady-state universal 
relation in Fig|| and finally comes back to the dashed 
curve. The deviation from the steady-state curve be- 
comes larger with increasing 7, as in the experiments of 
elongational flow. 



VI. TIME-CORRELATION FUNCTIONS AND 
TUMBLING IN SHEAR FLOW 

In Fig|ll], we show the end-to-end vector correlation 
functions C(t) = (P(t) ■ P(0))/(|P| 2 } with and with- 
out shear flow for various 7 at T — 0.2. For 7 ^ 0, 
it rapidly decreases, negatively overshoots, and finally 
approaches zero with dumped oscillation superimposed. 
This oscillatory behavior arises from random rotation of 
chains in shea*, flaw, which is well known in dilute poly- 
mer solutionsc3c3 but has not been reported in polymer 
melts. This is more evidently seen in Figjl^, where we 
show time development of the x component of the end- 
to-end vector Pj = Rn — Ri of one chain for 7 = 10~ 3 
(a), 10" 2 (b), and 10" 1 (c) at T = 0.2. The correspond- 
ing Weisenberg number (4.1) is given by Wi = 60,600, 
and 6000, respectively. In Fig.|l2](d), we show chain con- 
tours projected onto the x-y plane at the points 1 ~ 8 
indicated in Fig.[l2](c). When the chains change their ori- 
entation, their shapes are contracted as in the case of a 
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single chain in solution£3 The average period of tumbling 
is about 35/7 in our case. 

We may introduce the van Hove time-correlation func- 
tion (2.14) even in shear flow if the particle displacement 
vector is redefined asEl 

ARj(t) = Rj(t + t ) - Rj(t ) 

-7 / dt'Y G (t +t')e x , (6.1) 
Jo 

where Yq is the y-component of Rq = N~ lJ £2^ =1 Rj. 
From the net displacement, the first term, we have sub- 
tracted the flow- induced displacement, the second term. 
Figure |l^ shows F q (t) with q = 2tt for various 7 at 
T = 0.2. Comparison of this figure with Fig.|| suggests 
that applying shear is analogous to raising the tempera- 
ture. This tendency was already, reported for the case of 
supercooled binary mixtures oE-j 

We introduce the shear-dependent Rouse time tr = 
tr(T, 7) and the a relaxation time T a = r Q (T, 7) by 

C(r R ) = e- 1 , F g (r ce )=e- 1 . (6.2) 

We may then examine how shear can accelerate the mo- 
tions of chains and individual beads in shear flow. Fig- 
ure [l4| shows tr and r a as functions of 7 at T = 0.2. 
In our short chain system, both tr and r a decrease for 
7 £ ^{T^O)- 1 < TaCT.O)- 1 . Our data at T = 0.2 are 
consistent with 

mC^T)" 1 = m (T,o)- 1 [i + A R7 ] ) (6.3) 

T a (T, 7) _1 = Ta (T, 0)- 1 !! + (A a j)% (6-4) 

where A R 10 4 ~ t r (T,0), A a = 6000 ~ 20t q (T,0), 
and \x = 0.77. The average tumbling period in Fig. 12 is 
about 4tr(T, 7). For simple supercooled liquids we al- 
ready introduced the van Hove time-correlation function 
in shear flowO and obtained r Q (T,7) _1 = T a (T, 0) 1 [1 + 
A a7 ] with A Q ~r a (T,0). 

The sensitive shear-dependence of r Q (T, 7) predicted 
by (6.4) suggests potential importance of dielectric mea- 
surements in shear flow.cl As a first experiment, Mat- 
suyama et al. measured the dielectric loss function 
e"(ti;,7) in steady shear 7 in oligostyrene and polyiso- 
prene meltsJi3 In the former melt at T = 42°C, e"(w,7) 
decreased nonlinearly as a function of 7 at low frequencies 
(u; < 10 5 s _1 ) in the shear-thinning regime (7 > 15 s" 1 ). 
Their finding indicates that T a decreases as a function of 
7 in the non-Newtonian regime, consistently with (6.4), 
More systematic dielectric measurements in supercooled 
systems under shear flow are very informative. 

VII. SUMMARY 

We have performed very long MD simulations of a su- 
percooled polymer melt composed of M = 100 short 



chains with bead number N — 10 in quiescent and 
sheared conditions. We here summarize our main sim- 
ulation results together with remarks. 

(i) The stress relaxation function G(t) is shown to fol- 
low a stretched exponential decay (1.1) on the scale of 
the a relaxation time r Q and then the Rouse relaxation 
(3.6) on the scale of tr. 

(ii) The nonlinear shear regime sets in at extremely 
small shear rate of order t^ 1 in supercooled states, where 
marked shear-thinning and shape changes of chains are 
found. Scattering and birefringence experiments from 
weakly sheared melts near T g seem to be very promising. 

(iii) In the nonlinear shear regime, each chain under- 
goes random tumbling in our melt as in the case of iso- 
lated polymer chains in shear flow. It is of great interest 
how this effect is universal in solutions and melts and how 
it influences macroscopic rheological properties. For ex- 
ample, we are interested in whether or not such tumbling 
occurs in sheared entangled polymers. 

(iv) Transient stress divided by 7 after application of 
shear flow obeys the linear growth J Q dt'G(t') ~ Got 
for strain less than 0.1 and then saturates into a non- 
Newtonian steady-state viscosity. This initial growth is 
much steeper than that predicted by the Rouse model. 
As a result, the stress-optical relation does not hold tran- 
siently under deformation even in the linear (zero-shear) 
limit. Its violation is more enhanced for larger shear 
rates. These are consistent with the experiments. 

(v) The time-correlation functions in shear flow are cal- 
culated for the end-to-end vector and the modified par- 
ticle displacement in (6.1). The former represents the 
relaxation of chain conformations, while the latter the 
monomeric relaxation on the spatial scale of the particle 
distance. We can then determine the shear-dependent 
relaxation times, tr(T, 7) and T a (T, 7). They decrease 
nonlinearly and behave differently as functions of 7 in 
the nonlinear shear regime as in (6.3) and (6.4). It is 
then of great interest how these times behave in strong 
shear for much larger N. We conjecture that if N is 
sufficiently large, shear should first influence the overall 
chain conformations, while it does not much affect the 
monomeric relaxations. We propose dielectric measure- 
ments in— shear flow, which should give information of 
t q (T, 7) £3 In addition, as demonstrated in Fig. 13, the 
effect of shear on the van Hove self-correlation function 
is analogous to raising the temperature above T g as in 
supercooled binary mixtures J§EJ 
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( 

FIG. 1. Normalized end-to-end vector time- correlation 
function C(t) in (3.4) for T = 1.0, 0.4, and 0.2. The dot- 
ted lines are the results of the Rouse model (3.5). The Rouse 
time tr in (3.2) is indicated by arrows. 




t 

FIG. 2. The van Hove self-correlation function F q (t) at 
q — 2tt for T = 1.0, 0.4, and 0.2 on a semi-logarithmic scale. 
The dotted line represents the stretched exponential decay 
oce X p[-(t/r a ) - 64 ]. 
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t 

FIG. 3. The stress relaxation function G(t) in (3.10) 
(thin-solid lines) at T — 0.2 in a supercooled state and T = 1 
in a normal liquid state. For T = 0.2, it can be fitted to 
the stretched exponential form, exp[— (t/r B ) "'] with r s = 90, 
(dotted line) for 1 < t < 10 3 and tends to the Rouse relax- 
ation function Gr(£) (bold-dashed lines) at later times. 




t 

FIG. 4. Comparison of the time-correlation function Gb{t) 
in (3.14) (thin-solid lines) and the Rouse relaxation function 
Gr(£) in (3.6) (bold-dashed lines). 
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FIG. 6. Universal stress-optical relation a xy /T vs Q xy in 
steady states under shear flow for T = 0.2, 0.4, and 1. 
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FIG. 7. (a) Isointensity curves of g s (r) in (4.3) in the x-y plane (—3.75 < x,y < 3.75, z = 0). (b) Those of the incoherent 
structure factor S(q) in (4.4) in the q x -q y plane (— tv < q x ,q y < tv, q z = 0). The values on the isolines are 0.01 + 0.02n in (a) 
and 0.1 + 0.2n in (b) with n = 0, 1, • • • ,4 from outer to inner. Here T = 0.2, 7 = 10~ 4 , and the flow is in the horizontal (a;) 
direction. The is the angle between the average chain shapes and the y axis. 
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FIG. 8. tan 8, 1 — 01/02, and Q xy vs 7r R in steady states. 
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FIG. 9. Shear stress divided by shear rate o- xy (t)/j vs t/TR 



for 7 = 10 _1 , 10" 



10" 



10 -4 (thin-solid lines) at T = 0.2 



where tr = 6 x 10 4 . The curves follow the linear viscosity 
growth function (bold-solid line) for -yt < 0.1, but depart 
from it for •yt > 0.1. The linear growth function in the Rouse 
model is also plotted(bold-dashed line). The arrows indicate 
onset of the nonlinear behavior. 
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FIG. 10. Parametric plots of a xy (t)/T vs Q xy (t) after ap- 
plication of shear at t — for T = 0.2. The curves initially de- 
viate from the universal steady-state curve obtained in FigJ^ 
(dashed line) but approach it ultimately. The deviations in- 
crease with increasing 7. 
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FIG. 11. Normalized time-correlation function of the 
end-to-end vector C(t) in (3.4) at T = 0.2 for 7 = 0, 10" 4 , 
10" 3 , 10" 2 , and 10" 1 from right to left on a semi-logarithmic 
scale. The negative overshoot for 7 > arises from rotational 
motions of chains. 
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FIG. 12. Time- evolution of the x component of the end-to-end vector P x (t) = Xn (t) — Xi (t) of one chain vs jt. Here T = 0.2 
and 7 = 10 _3 (a), 10 _2 (b), and 10 _1 (c) from above. Typical tumbling motions at the points 1 ~ 8 indicated in (c) are shown 
in (d), where the chain conformations are projected on the x-y plane. 
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FIG. 13. The van Hove self correlation function (3.7) with 



(6.1) at T = 0.2 for 7 = 0, 10 , 10 , 10" 
right to left on a semi- logarithmic scale. 
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FIG. 14. Two relaxation times tr(7) and r a (7) as func- 
tions of shear 7 at T = 0.2 determined from (6.2). Both these 
times decrease for 7 > tr(0) _1 ~ N~ 2 T a {0)~ 1 in our short 
chain system. The solid and dashed lines represent (6.3) and 
(6.4), respectively. The slopes of the curves at high shear are 
— 1 for tr(7) and —0.77 for 7-^(7). 
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